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Abstract. Geometric details of a nuclear reaction zone, at the time of par- 
ticle emission, can be restored from low relative- velocity particle-correlations, 
following imaging. Some of the source details get erased and are a potential 
cause of problems in the imaging, in the form of instabilities. These can be 
coped with by following the method of discretized optimization for the restored 
sources. So far it has been possible to produce 1-dimensional emission source 
images, corresponding to the reactions averaged over all possible spatial direc- 
tions. Currently, efforts are in progress to restore angular details. 
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1. Introduction 



In astronomy, interferometry was initially just used for determining star diameters 
or binary-star separations. Since then, the methodology has moved to the deter- 
mination of detailed images of star systems evolving with time, as exemplified in 
Fig. ^ This concerns the intensity and phase interferometry. It is needless to say 
that the nuclear reaction physics is way behind the astronomy, in extracting in- 
formation from particle interference or, generally, correlations. Besides astronomy, 
widely ranging areas, that employ different types of imaging, include police work, 
tomography and seismography. 

The general task of imaging can be formulated rather simply. Principally, it 
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CAP ELLA: COAST IR at 1290nm on 25/10/97 
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Fig. 1. Reconstructed image of Capella, from data taken on October 25, 1997, at 
wavelength of 1.3 /im, after 

amounts to inverting an integral relation of the form 



Here, C represents the quantity that is measured as a function of q and S, which is 
a function of r, represents an image of interest. Given the data on C, with errors, 
one needs to determine S, which requires an inversion of the kernel K. 

The possibility of imaging in the reactions of heavy nuclei, with many-particle 
final states, results from an interplay between the geometric features of the emission 
zone and of structures in two-particle wavefunction, usually found at low interparti- 
cle velocities, cf. Fig.|5J The features in the wavefunction normally change rapidly as 
the interparticle velocity changes. If we consider an inclusive cross-section for pair 
emission, da/dp 1 c?p2, it is given in terms of a multiparticle amplitude squared, inte- 
grated over unobserved particles. This amplitude can be factored into the outgoing 
two-particle wavefunction and a reminder. The reminder squared and integrated 
over unobserved particles represents the source S", without a strong dependence on 
the relative velocity or momentum q: 



Here, r and P are the relative cm. position and total momentum, respectively, and 
S' may be interpreted as a distribution of emitted pairs in r. 

If we normalize the two-particle cross section to the product of single-particle 
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Fig. 2. The possibility of imaging, in the reactions of heavy nuclei, results from 
an interplay of the coarse geometric features of the reaction zone, with pronounced 
and calculable geometric features within the two-particle wave-function. 

cross sections, we obtain the correlation function: 



in terms of the source function S now normalized to 1, J dr Sp(r) = 1, since, for 
large g, the correlation function will approach 1. For a plane- wave wave-function, 
we have l^l 2 = 1 and expect C ~ 1. On the other hand, any nonuniformities in 
|$| 2 can be used as a lens for testing S. 

The source function is subject to different expectations depending on conditions 
in a reaction. Thus a rapid freeze-out, with not much collective motion developed 
in a reaction, should lead to a compact source that is roughly spherical in shape. A 
rapid freeze-out, combined with a strong collective motion, may lead to an oblate 
shape perpendicular to the velocities of emitted particles. If many of the measured 
particles stem from secondary decays, the source function will exhibit an extended 
tail. If the particles are emitted from a long-lived residual system, the source 
function will develop a smutch shape along the particle velocities. By studying the 
source shape, one can learn about the nature of the reaction. 

2. Imaging for the Reactions 

In practice, information is only contained in the deviation of the correlation function 
from unity. If we subtract 1 from the sides of Eq. [31 we get an imaging relation 




da 




ft P (q) = Cp(q)-l= / dr(|^-)(r)| 2 -l) S P (r) = / drK(q,r) 5 P (r) . (4) 
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Source imaging is possible when |$q ^(r)| 2 deviates from 1, as due to symmetrization 
or interaction within the pair. 

Upon spin averaging, the kernel K depends only on the relative angle between q 
and r. This is important in the angular decomposition of S and 1Z and in the decom- 
position of the imaging relation. Upon decomposing the kernel and the functions 
according to 

K(q, r) = £)(2A + 1) K x (q, r) P x (cos 9) , (5) 

A 

and 

TZ{q) = V4^^ft Am (g)Y Am (q), S(r) = V^J^ ^"W Y Am (r) , (6) 

Am Am 

we get, from the relation Q, 

K P m (q) =4njdr r 2 K x (q, r) S Am (r) . (7) 

The above shows that different multipolarities of deformation for the source and 
correlation functions are directly related to each other. The A = version of the 
relation is 

K P (q)=4TT [ drr 2 K (q,r)S P (r), (8) 



where IZ(q), K and S(r) are the angle-averaged correlation, kernel and source, 
respectively. 

The simplest case of imaging corresponds to a situation where the structure of 
|$| 2 is due to a pure interference, e.g. neutral pions or gammas. For a symmetrized 
wavefunction, 

$<-)(r) = i=(e^ r + e - 4qr ) (9) 

the kernel is 

K(q,r) = |$£->(r)| 2 - 1 = cos (2q • r) , (10) 
and the imaging amounts then to a Fourier-cosine transform 

TMq) = J drcos(2qr)Sp(r) =*> S P (r) = i J dq cos (2qr) K P (q) . (11) 

The angle-averaged relation also represents a Fourier transformation, since 

sin (2qr) 

^ 0(<Z ' r) = ^^' (12) 

and, in general, the relations for different multipolarities represent Fourier-Bessel 
transforms, since K\(q,r) — (— 1) A / 2 j\(2qr). As an example of the application 
of Fourier transformation in imaging, Fig. [3] shows the relative source of negative 
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Fig. 3. Relative source of negative pions in Au + Au collisions at 10.8 GeV/nucleon, 
from Fourier transformation of the negative-pion correlation function obtained by 
the E877 Collaboration [Q. 

pions obtained from the Coulomb corrected correlation function from measurements 
of the E877 Collaboration 

For many particle pairs, as e.g. pp, the interaction can be neither ignored nor 
simply corrected for, and the straightforward Fourier transformation cannot be 
used. In the general case of correlation, though, the source discretization with a x 2 
fitting may work. For this, e.g. in the angle-averaged relation, one discretizes the 
integral to get 



Hi = J2 An Ar r2 J K o(Qi,rj) S( rj ) = Kij Sj , 



(13) 



and, subsequently, varies the Sj values to minimize 



E 



(14) 



The derivative of \ 2 with respect to Sj gives a set of algebraic eqs. for S: 

Y / ^(K ij Sj-Tir P )K i j = 0, 



(15) 



with the solution, in a matrix form, 

S = (K T K)- 1 K T n exp . 



(16) 
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Figure ^ shows the results of simple tests concerning imaging, where a cor- 
relation function is constructed out of an assumed source, errors are added and 
inversions are carried out for the source. This figure illustrates the problems that 
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Fig. 4. On the left, pp correlation functions constructed from two assumed source 
functions, without (lines) and with (symbols) errors added. On the right, the re- 
sults of source restoration through imaging, from the sharper of the two simulated 
correlation functions on the left, following either the discretization with an assumed 
bin size of 2 fm (right top figure) or the optimized discretization method [0 (right 
bottom). 

can be generally encountered in imaging. Specifically, the correlation function may 
not reflect fine details of the source. If excessive detail is demanded within some 
region, the inversion of the kernel begins to remind the inversion of a zero, which 
is signalled by large changes in the output for small changes in the input. This 
situation is encountered in the naive imaging in Fig. ^ relying on a fixed 2 fm bin 
discretization of the seeked source. The situation may be remedied by following the 
method of optimized discretization, where the source discretization is varied prior 
to the inversion, to ensure minimal anticipated errors. An exemplary result of the 
procedure, together with the original source, is also further shown in Fig. ^ 
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The imaging has significantly changed the interpretation of different pp corre- 
lation data previously described in terms of Gaussian sources [0], cf. Fig. [3] Thus, 
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Fig. 5. On the left, pp correlation functions measured for different indicated total 
momentum ranges, in the 14 N + 197 Au reaction at 75 MeV/nucleon. On the right, 
the sources imaged (shaded regions) from the correlation functions on the left, 
together with the results (lines) when assuming sources in the Gaussian form, after 



the short-range portion of sources in the reaction is, generally, non-Gaussian and 
the sources exhibit long-range tails. The weight of the preequilibrium short-range 
portion of a source, relative to the long-range halo, changes strongly with the mo- 
mentum of an emitted pair, while the spatial extent of the preequilibrium portion 
remains relatively constant with the momentum. The changes in the relative weight 
of preequilibrium particles get reflected in the height of the measured correlation 
function. The changes in the height of the correlation function get, though, mis- 
interpreted, in terms of the changes in extent of the preequilibrium shape, when 
forcing the Gaussian shape onto a source. 



3. Towards 3D Imaging 

Anisotropics in correlation functions have been used to access differences in typical 
emission times for different particles [ 5 and to access asymmetries in the emission 
zone from reaction geometry and collective motion [ |Uj . The basic relation that 
can be used to access the anisotropics of the source within imaging is Eq. J7J . For 
identical particles, only even A contribute to the angular decompositions in (|f f (I 
and Unj). 
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For the decomposition, we need to choose a reference coordinate system where 
asymmetries are assessed and a convenient system is one with the z-axis along the 
beam axis and the x-axis along the net transverse momentum of the emitted pair. 
Either the source S or the correlation C (or TV) may be represented in that system 
as 

C(q) = E/ ^'C(q')VX m (q')n m (q) = ^ f dtf C(q') Re (Y xo (q')) 

Am A 

xRc(y A0 (q)) + 2 [ dQ.'C(q')Re(Y Xm (q')) Rc(r Am (q)) 

= E E ^(g)^-..^ = / ^ W) + 3 / f£ C«) 

A fci . . .k\ 

where the symmetry with respect to the y — > —y reflection and the real nature of 
C have been exploited after the second equality. (The additional symmetry in a 
symmetric system at midrapidity additionally limits the contributions to only terms 
characterized by even A + to.) The expression after the third equality is in terms 
of cartesian components for spherical tensors of rank A. Finally, the last expression 
makes the expansion in terms of tensors of increasing rank explicit. 

The leading term on the r.h.s. of l|17|l is the correlation function averaged over 
angles, C^°\q). The dipole tensor in the second term on the r.h.s. of (|17|l . i.e. a 
vector, may be expressed as 

7f (9)=3| ^C(q')^ = C«( g )n«(g); (18) 

the dipole deformation is characterized by a direction angle 9^ and deformation 
parameter , both dependent on q. The quadrupole tensor in the third term can 
be next expressed as 

= Cf W^ [2) +C^WW {2) -{C? + C?)m (2) . (19) 

The quadrupole deformation may be characterized in terms of two deformation 
parameters, such as the larger C^ 1 and lower c[ 2 \ within the P^-beam plane, 
and in terms of the direction angle for the larger deformation, 9^ 2 \ The axis 2 for 
the deformation tensor, by symmetry, is perpendicular to the P^-beam plane. An 
analogous decomposition holds for the source. Overall, for deformations in terms of 
the lowest angular moments, there are 6 parameters, as a function of q or r, needed 
in describing the correlations or sources for nonidentical particles, and 4 parameters 
for identical particles, since the dipole deformations then vanish. 
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In terms of numerical analysis, for A > 1, the same methods can be employed 
as for A = 0, including basis splines and the optimized discretization. The novel 
issue encountered for A > 1 is the theoretical deterioration of angular resolution as 
q — ► or r — ► oo, amplified by the possible experimental deterioration of resolution 
for any practical data. 

Figure 6 shows results from a test similar to that in Fig. 4, where now a 3- 
dimensional correlation function has been constructed from an assumed source for 
7r~ pairs The function has been modified by errors, an angular decomposition 
was applied, and then inversion was employed for each multipolarity and the source 
was reconstructed. The original and reconstructed sources are compared along 
three orthogonal directions in Fig. 6 and it is seen that the imaging is successful at 
a semiquantitative level. 




Fig. 6. Comparison of the imaged (bands) and input (symbols with errors repre- 
senting finite event statistics) 7r~ — 7r~ sources along three perpendicular directions, 
after [0. 



4. Conclusions 

We have demonstrated the feasibility of imaging the emission sources in heavy-ion 
reactions. The one-dimensional imaging has been carried out for data on pion, pro- 
ton, kaon and IMF correlations. The three-dimensional imaging is under practical 
developments and we have here described its essential elements. The methods that 
are advanced, in particular involving optimized discretization, permit to investigate 
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source images, on a logarithmic scale, for large separations. The imaging gives 
access to the emission source structure. In their details, the images can contain in- 
formation pertaining to the phase-space density of particles at freeze-out, produced 
entropy, spatial freeze-out density, reaction geometry, developed collective motion, 
resonance decays and prolonged emission including phase-transition effects. 
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